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Abstract. - We use an off - lattice bead - spring model of a self - avoiding polymer chain 
immersed in a 3-dimensional quenched random medium to study chain dynamics by means of 
a Monte - Carlo (MC) simulation. 

The chain center of mass mean-squared displacement as a function of time reveals two 
crossovers which depend both on chain length N and on the degree of Gaussian disorder 
A. The first one from normal to anomalous diffusion regime is found at short time Ti and 
observed to vanish rapidly as t\ oc A~ n with growing disorder. The second crossover back 
to normal diffusion, ti, scales as T2 oc N 2lJ+1 f (N 2 ~' iv A) with / being some scaling function. 
The diffusion coefficient Dn depends strongly on disorder and drops dramatically at a critical 
dispersion A c oc ]\f- 2 + Sv Q f the disorder potential so that for A > A c the chain center of mass 
is practically frozen. 

The time-dependent Rouse modes correlation function C p (t) reveals a characteristic plateau 
at A > A c which is the hallmark of a non - ergodic regime. These findings agree well with our 
recent theoretical predictions. 



Introduction. Considerable progress has been made in understanding of the static 
[TlEll5llll31lrJir71l51IUl [Tni[n] and dynamic [TMHflin] behavior of polymer chains in a quenched 
random medium. It was shown already in the paper by Cates and Ball |2 that for sufficiently 
large diffusion coefficient (so as the chain may sample the whole system) the quenched and 
annealed problems are equivalent. Indeed, since the chain always moves towards the deepest 
energy minimum, quenched and annealed problems can differ only for finite volume of the 
random medium |10U11| . 

In a recent publication [15) we have argued that dynamic aspects can significantly affect 
the problem in question. We have treated this problem |15| using the generating functional 
method which based on the Langevin dynamics. Within this method the averaging over 
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the quenched disorder leads to coupling between different dynamic trajectories which in turn 
eventually results in a non - Markovian diffusional slowing down, i.e. in anomalous diffusion, 
as well as in non - ergodic behavior for the Rouse modes of the polymer chain. The resulting 
equation of motion (which was derived within the Hartree approximation) has a memory 
kernel typical for the mode - coupling theory of the glass transition ^j- Self - consistent 
treatment of this equation leads to the conclusion that in d - dimensional space the center of 
mass diffusion coefficient decreases according to the law: 

D N w D R [1 - const(A/6 d )iV Q ] , (1) 

where A is the dispersion of the disorder potential, b is the Kuhn segment length, and a = 
2 — vd (y is the Flory exponent). This differs remarkably from the earlier prediction, Dn ~ 
Z?r exp[— (A/b d )N a ], based on a simple Markovian diffusion argumentation O0J^]. As a 
result some critical disorder, A c « b d N~ a , exist such that at A > A c the chain center of mass 
is arrested by the random potential. In the case when the capture time i cap >> i oq (where t eq 
is the time of internal degrees of freedom equilibration) the quenched and annealed problems 
are not equivalent even for an infinite volume! Rouse modes of the chain also freeze at a 
common disorder strength A, so that the non - ergodicity plateau shows up according the A - 
type (or continuous) transition JHj. In the present letter we will mainly check these important 
predictions of the ref. [15] by the direct and vast MC - simulations. In the conclusion we also 
sketch some possible future developments. 

Model. - In the present study we employ an off-lattice Monte Carlo model of a bead- 
spring polymer chain, used in numerous investigations before ^Jj. During the simulation 
beads are chosen at random and displaced to a new trial position in their environment subject 
to a Metropolis procedure. Each polymer chain consists of N beads, connected by anharmonic 
springs which are described by the finitely-extensible nonlinear elastic (FENE) potential JS| 

Ufene = ^ — ^ n 

where I is the bond length which can vary in between l m - m < I < ^ max , and has the equilibrium 
value ioj with R = Z max — I, l m in = 21 — imax> and K is the spring constant. As before |18j . 
we take ^ max = 1 as the unit length, / m ; n = 0.4, Iq — 0.7, and K/ksT = 40. The non-bonded 
interactions between the beads are described by the Morse potential ^7] 

Um{t) = em exp [-2a (r - r min )] - exp [-a (r - r min )] , (3) 

where r is the distance between the beads, r m i n = 0.8, cm = 1> and a = 24 (since the potential 
has a very short range, the simulation code is rather fast ^E])- The 0-temperature for the 
model is ksTe ~ 0.62, so that at ksT = 1, used throughout in this study, the chain is in a 
good solvent regime and all energetic variables are dimensionless. 

The energetic disorder is created by dividing the simulation box into 64 3 hypothetical 
subcells, each of size l max = 1, which are ascribed a Gaussian-distributed local potential so 
that its dispersion is taken as a measure of the degree of disorder A. Typically runs of length 
10 s -i- 10 9 MC steps (MCS) have been carried out and observables have been averaged over 
500 different realizations of disorder. 

Results. - First of all we emphasize that the stochastic dynamics in the system under 
study is not time translational invariant (TTI) |19|. This means that in the time - dependent 
correlation function one should keep control on two times: the waiting time t w , i.e. the time 
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elapsed after the system preparation (aging) before any measurement starts, and the actual 
measurement time t. With this in mind we have studied the center of mass mean - square 
displacement (MSQD) 

g 3 (t, t w ) = ([r CM (t) - r CM (iw)] 2 ) , (4) 

where the center of mass position tcmW = J2iLi and the angular brackets denote 

an average over the realizations of disorder. In contrast, in cases of vanishing disorder TTI 
holds and the time correlation functions depend on a single variable St = t — t w . 
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Fig. 1 - Mean - square displacement of a central monomer, gi(t), of an end monomer, g&(t), and of 
the chain center of mass, gz(i) vs. time for N = 32 and A = 0.95. Dashed lines with slope 1 denote 
normal diffusive behavior whereas the crossover times ri (from) , and T2 (back to) normal diffusion 
are marked by arrows. The mean squared radius of gyration of the polymer, R s (t), appears as a 
horizontal line. 

PigHJshows a typical example for the chain length N = 32 the disorder strength A = 0.95 
and the MC - time interval t > 10 8 . One can clearly distinguish three diffusional regimes. In 
the initial short interval < t < t\ a relatively fast normal diffusion shows up. The crossover 
time n corresponds to the moment when disorder starts to effect the chain motion. In ref. |15) 
this moment was denoted by and was defined as the point where the memory - friction term 
overwhelms the bare - friction term. After that the non - Markovian diffusion gets under way. 
In the interval t\ < t < T2 the MSQD manifests itself as anomalous diffusion with gz{t) oc t@, 
where the exponent (3 < 1. At time t > T2 (which corresponds to the MSQD g$ (t) ^S> Rg) the 
diffusion goes back to normal albeit with a significant slowing down (see below). 

The MSQD - plot given in FigQJhas been obtained for waiting time t w = 0. In the present 
report we shall not consider aging phenomena (i.e. the t w - dependence of the time correlation 
functions as well as the violation of the fluctuation - dissipation theorem (FDT)) putting this 
off for a later publication. Nevertheless FigEl clearly demonstrates the drift of the full energy 
E with time t which for N = 8 gets pronounced at A > 0.9. Evidently in this case the TTI 
is violated for larger disorder. The age of the system also influences static conformational 
properties as the mean - squared radius of gyration Hz of the polymer chain. Figure [3] shows 
that at small disorder, A = 0.25, the i?g versus TV dependence has the same slope as in the 
pure (i.e. without any disorder) system. An inspection of the measured Rl vs. N relationship 
yields v = 0.605 ± 0.005 for the Flory exponent. At large disorder, A = 1.00, however, the 
slope of R g (N) is affected by the waiting time t w . The value of t w which we have used in 
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Fig. 2 - Energy (per monomer) vs time of a chain with N = 8 for growing degree of disorder A (given 
as a parameter). 

Fig. 3 - A scaling plot of the radius of gyration, E? g , with chain length N for small, A = 0.25, and 
large, A = 1.00, degree of disorder. 



FigED corresponds to the slow diffusion regime, i.e. t w > T2- In this case v = 0.68 ± 0.01, 
i.e. v > va=o which is consistent with the results obtained on hierarchical lattices for the 
quenched problem. The quenched problem in ref. corresponds to the case in which the one 
end of the chain is kept fixed. In this case the polymer stretches away from the fixed end in 
order to find a region with a minimal energy. The Imry - Ma arguments have been used in 
ref [Sj to show that f qU cnchcd = 2/3. 

Now we go back to MSQD - plot (Fig|TJ to discuss how the crossover times, t\ and ra, 
depend on A. This is presented in Fig0] 
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Fig. 4 - Crossover times t\ and T2 (inset) for chains of length N at different degree of disorder A. 



The first crossover n decreases whereas the second one t 2 increases with the disorder 
strength A, so that the anomalous diffusion interval is extended. While for the relatively 
long chains, N > 32, the initial interval n is getting progressively so short that a scaling 



A. MlLCHEV, V.G. ROSTIASHVILI and T.A. VlLGIS : DYNAMICS OF A POLYMER IN A QUENCHED RANDOM MEDIUM: A MONTE CARL 



dependence on A is difficult to establish, for shorter chains (where the statistics is better too) 
the Ti vs A relationship can be reasonably fitted by t% oc A -11 . In ref. JTjjjJ a theoretical 
prediction for n was derived as: 

nocA-W-fl , (5) 



where the anomalous diffusion exponent /? = 1 — (2— vd)/(2v + 1) < 1. In the case of d = 3 
and = 0.6 one obtains j3 = 0.9 which leads to T\ ~ A -10 . This prediction agrees closely 
with the above - mentioned MC findings. As is clearly seen in the inset of Fig. all data 
for the second crossover time T2 follow a scaling law T2 oc N 1+2v f(AN 2 ~ 3l/ ) where / is some 
scaling function. 

As mentioned above, at t > r 2 the normal diffusion regime shows a significant slowing 
down. This extremely important result is depicted in FigJSJover a wide range of chain lengths. 




Fig. 5 - Variation of the normalized diffusion coefficient Dn/Dr with scaled degree of disorder 
x — AiV 2 ~ 3l/ for chains of length 8 < N < 256. Dashed line denotes the simple Markovian diffusion 
prediction NDn oc exp(— x) 4 15 , a full line is a best fit with NDn oc (x c — x) ' 8 . Unsealed data is 
shown in the inset. 



All points in the reduced coordinates (NDn vs AN 2 ~ 3u ) lie on a single master curve. It ap- 
pears that our MC - data nicely follows the best fit Dn oc N~ 1 (x c — x) ' 8 (where x = AN 2 ~ 3l/ ) 
which correlates well with the theoretical predictions for Dn (see eq.Q). A remarkable and 
strong systematic difference between the master curve (solid line) and the simple Markovian 
diffusion result (dashed line) can be seen, which grows progressively with increasing of degree 
of disorder. The basis for this slowing down is the coupling between diffusional trajectories 
(or between different parts of the same trajectory) caused by the averaging over the disorder. 
This furnishes the non - Markovian slowing down regime, which mathematically shows up as 
a memory kernel term in the resulting equation of motion |15) . Such behavior is generic for 
the "schematic" mode - coupling theory of the glass transition |16|. The existence of such 
dynamic singularity for a finite chain is in principle questionable. More precisely one should 
better mean that the capture time of the center of mass, t cap , is finite albeit astronomical, so 
that in the slowing down regime t cap ^> exp(AiV 2-3l/ ). This time scale goes beyond the 
scope of our present MC - simulation. 

So far, we have focused mainly on the chain center of mass slow dynamics. We studied 
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also the time dependent correlation function of the Rouse modes [20] 

Cp(Mw) = {X p (t)X p (t w )) , (6) 

where the Rouse mode component with an index p, X p (t) = N^ 1 Yli=i r i(t) cos[(i— 0, 5)p7r./V -1 ],| 
is measured at time t. This correlation function reflects the dynamics of the internal degrees 
of freedom. It can be seen in Fig. [HJthat at relatively large disorder strength A a characteristic 
plateau appears which is reminiscent of a non - ergodic behavior for the so called A - type 
dynamical transition J^j. Figure El illustrates this non - ergodic plateau behavior for the first 
mode (p = 1) and the chain length TV = 16. 
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Fig. 6 - Variation of the plateau height h with degree of disorder A of the Rouse-mode correlation 
function Ci(t,0) of a chain with N = 16. In the inset Ci(t,0) are shown for the respective values of 
A. 

The hight h of the plateau is found to be a continuous function of A, i.e. h oc (A— 1.06) 2 ' 5 , 
which qualitatively agrees with the recent theoretical predictions |15) . It is interesting to note 
that the critical A c w 1.06 coincides with the critical value A c which is obtained for N = 16 
from the diffusion results (see the FigJSJ. The well - known F± model, which has been studied 
in ref. |16j . provides a simple example of such kind of continuous transition. 

Conclusion. - In this work an off - lattice model of a self - avoiding polymer chain in 
quenched random field is used to study the slow polymer dynamics. We examine the center 
of mass MSQD at different degree of disorder and find a broad regime of anomalous diffusion 
which starts at a crossover time t\ and ends at a second crossover T2 by going into normal 
diffusion regime. We have demonstrated that these crossover times are observed to scale in 
agreement with the theoretical predictions |15j . 

The diffusion coefficient Z?at(A) at t > drops dramatically at a critical value A c oc 
Simulation data for Dm of different polymer chain length and different degree of 
disorder is observed to collapse on a single master curve as predicted in a recent theory [To] . 
This significant slowing down effectively originates from the coupling between different dy- 
namic trajectories of a polymer chain caused by the averaging over many disorder realizations. 
The Hartree approximation which has been used in JHj is close in spirit to the well- known 
mode - coupling approximation |16| and provides the basis for the treatment of non - Marko- 
vian regimes. In contrast, a treatment in terms of Markovian arguments yields an exponential 
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relationship between diffusion coefficient and reduced disorder which proves to be inadequate. 
The same reasons bring about the Rouse modes freezing at a common disorder strength A c , 
whereby the generic non - ergodicity plateau appears continuously. This is also in line with 
the theoretical predictions 

Finally we should stress that the system at issue does not satisfy TTI so that the waiting 
time t w is a very important parameter. Although t w has been controlled in the calculation of 
MSQD (eq.Ql) and the Rouse mode correlators, eq.@), its impact on the dynamics is to be 
reported in a following communication. To this end the generating functional approach, which 
has been used in ref. |15|. should be generalized to the out - of - equilibrium dynamics |19) 
where the influence of i w (or the aging) and FDT violation are the effects of a primary interest. 
We plan to consider these problems in a future publication. 
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